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Summary 


It  was  the  objective  of  this  program  to  determine  the  time-dependent  (viscoelastic) 
response  of  an  elastomer  at  various  temperatures  so  as  to  be  able  to  construct  the  uniform 
master  relaxation  response,  from  which  the  mechanical  responses  under  a  wide  range  of 
prescribed  load  or  defonnation  histories  may  be  deduced  at  any  reference  temperature. 

This  report  summarizes  work  accomplished  during  the  referenced  grant.  Early  efforts  were 
expended  on  re-establishing  test  facilities,  verifying  proper  operations  and  their  precision 
limits.  Particular  attention  was  devoted  to  the  precision  with  which  the  loads  and 
deformations  could  be  determined.  Strain  measurements  were  taken  for  monitoring 
purposes  by  means  of  the  digital  image  correlation  method.  Because  of  the  desire  to 
determine  also  the  time  dependent  Poisson  behavior,  digital  image  correlation  was 
attempted  in  both  the  axial  and  transverse  directions.  While  consistent  data  were  obtained 
for  strain  monitoring  in  the  axial  direction,  the  precision  with  which  lateral  strains  could  be 
determined  overlong  periods  of  time  was  not  sufficient  to  specify  the  volumetric  or  bulk 
behavior  in  the  rubbery  domain. 

The  report  is  divided  into  several  parts,  namely  a  summary  of  linear  viscoelasticity 
principles,  the  delineation  of  the  experimental  requirements  and  verification  of  the  test 
method  precision.  In  this  context  it  is  of  interest  to  note  that  test  machine  precision  was 
found  to  be  insufficient  to  prescribe  axial  strains  on  the  uniaxial  specimens  of  below  about 
1-2%.  However,  a  2%-strain  prescription  was  still  able  to  render  the  essentially  linearly 
viscoelastic  behavior  of  the  material. 

The  time-temperature  superposition  method  provided  a  master  relaxation  function  for  the 
elastomer  and  determined  its  glass  transition  temperature  to  be  near  -50°C.  The  range  of 
modulus  values  measured  over  a  temperature-reduced  time  range  of  15  decades  varied  by  a 
factor  of  about  12-15.  By  usual  standards  of  polymer  behavior  in  the  transition  range  this  is 
not  a  very  large  variation,  though  not  abnormally  small,  especially  if  one  considers  that  the 
full  range  of  the  glass-to-rubber  transition  has  not  been  captured.  Estimates  are  presented 
that  illustrate  the  conversion  of  the  measured  relaxation  response  of  this  material  into  a 
different  time  responses  under  short  term  loading  at  different  temperatures. 
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1  BACKGROUND 


1.1  LINEAR  VISCOELASTICITY 

The  present  investigation  was  limited  to  linearly  viscoelastic  behavior,  because  (a)  the 
understanding  of  non-linear  material  characterization  of  polymer  mechanical  behavior  is  in 
a  poorly  developed  state  and  needs  much  development  yet,  and  (b)  there  were  insufficient 
resources  and  time  to  embark  on  a  drawn-out  characterization  with  the  aim  to  develop  non- 
linearly  viscoelastic  material  characterization.  While  there  are  proposed  frame  works  for 
representing  nonlinear  material  behavior  (Schapery  1968,  1969),  Emri  and  Knauss  (1981, 
1987),  they  are  not  in  a  state  of  uniformally  accepted  engineering  practice.  The  framework 
to  describe  linearly  viscoelastic  material  behavior  used  effectively  for  many  years  in 
engineering  applications  bears  phenomenological  character.  It  is  based  mathematically  on 
either  an  integral  or  differential  formulation  with  the  material  representation  described 
realistically  in  numerical  (tabular)  or  functional  form(s).  The  fundamental  equations 
governing  the  linearized  theory  of  viscoelasticity  are  the  same  as  those  for  the  linearized 
theory  of  elasticity,  except  that  the  generalized  Hooke’s  law  of  elasticity  is  replaced  by  a 
constitutive  description  that  is  sensitive  to  the  material’s  (past)  history  of  loading  or 
deformation.  It  will  be  the  purpose  of  the  immediately  subsequent  sections  to  delineate  the 
details  of  this  formalism  of  material  description. 


1.1.1  Specific  Constitutive  Responses  (Isotropic  Solids) 

With  the  subsequent  and  specific  study  in  mind  the  following  discussion  is  limited  to 
isotropic  materials.  Recalling  that  the  stress  and  strain  states  may  be  decomposed  into 
shear  and  dilatational  contributions  (deviatoric  and  dilatational  components),  we  deal  first 
with  the  shear  response  for  exemplary  purposes  followed  by  the  volumetric  part.  Thermal 
and  pressure  sensitive  characterization  will  be  dealt  with  at  a  later  time. 


1.1.2  Shear  Response: 

Let  x  denote  any  shear  stress  component  and  s1  its  corresponding  shear  strain.  Denote  the 
material  characteristic  for  unit  step  excitation  in  the  present  shear  context  by  p(t).  This 
function  is  identified  as  the  relaxation  modulus  in  shear2  (for  an  isotropic  material).  The 
constitutive  response  in  a  pure  shear  deformation  is  the  given  by 

r(t)  =  2j//(t-£)^df  (1) 

o 


1  Standard  continuum  mechanics  notation  is  consistently  implied  in  this  document. 

2  Relaxation  is  a  phenomenon  associated  with  prescribing  a  step-wise  deformation  on  a  test  sample.  If  the 
deformation  mode  is  a  shear  strain,  the  resulting  relaxation  modulus  is  called  a  shear  relaxation  modulus.  It  is 
defined  as  the  resulting  time  varying  shear  stress,  normalized  by  the  amplitude  of  the  (step)  strain. 
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(2) 


=  2s(0)/u{i)  +  2  | //(t  - 

0J+ 

=  2/d  *  da 


where  the  star-notation  replaces  the  integral  notation  for  purposes  of  brevity  and  ju(t)  is 
clearly  the  relaxation  modulus  under  unit  step  strain  input  since  it  follows  from  (2)  that  for 
the  step  strain  history  a{t )  =  s0h(t)  of  amplitude  s0,  the  response  is 


r(t)  =  2p(t)  80 


(3) 


The  complementary  (inverse)  relation  for  a  prescribed  shear  stress  is 

=  (4) 

20  d £ 


2  2  Q+ 
1 

=  —J*cfr 
2 


(5) 


where  now  the  function  J(t)  is  the  shear  creep3  compliance,  which  represents  the  creep 
response  of  the  material  in  shear  under  application  of  a  step  shear  stress  of  unit  magnitude. 


1.1.3  Bulk  or  Dilatation  Response. 


Let  Sii(t)  represent  the  first  strain  invariant  and  iy(/)  the  corresponding  stress  invariant.  The 
latter  is  recognized  as  three  times  the  pressure  P(f),  i.e.  tu(/)=3P(/).  In  completely 
analogous  fashion  to  equations  1  and  2  the  bulk  behavior,  governed  by  the  bulk  relaxation 
modulus  k(t)  is  represented  by 

t  da--  ( (ft 

a jj  (t)  =  3 J  k(t  -  £)  — d^  (6) 

o  d^ 


=  3*ii(0)k(t)  +  3jk(t-£) 


o+ 


dgjj(^)  dg 

d£ 


(7) 


=  3k  *  ds{i 


The  inverse  relation  is,  similarly, 


(8) 


3  Creep  is  the  time  dependent  strain  response  of  a  test  body  to  the  prescription  of  a  time-invariant  stress.  The 
creep  function  is  a  monotonically  increasing  function  of  time. 
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(9) 


1  1  r  da ::({;) 

=  -ajj(0)M(t)  +  -  \M(t -S)  JJ  ^  d$ 

0+  ^ 

=  —  M  *dsij 
3 

where  the  function  M( t)  represents  now  the  dilatational  creep  compliance  (also  called 
the  bulk  creep  compliance),  or,  in  physical  terms,  the  time  dependent  fractional  volume 
change  resulting  from  the  imposition  of  a  unit  step  pressure. 


1.1.4  Uniaxial  behavior. 

It  is  very  common  in  engineering  descriptions  to  work  with  the  uniaxial  characterization 
commensurate  with  the  use  of  Young’s  modulus  in  linear  elasticity.  If  one  defines  E(t)  as 
the  relaxation  modulus  characterizing  the  uniaxial  stress  response  to  a  step  strain  input  of 
unit  magnitude,  then  the  stress  response,  say  on(t)  following  a  uniaxial  step  strain  history 
Sn(t)  is  gicen,in  complete  analogy  to  equations  (1-2)  as 

g~n  (0  =  f  E(t  -  g)  —~r^~  (10) 

o  ^ 


=  sn(0)E(t)+\E(t-4) 

o+  “S 

=  E  *  ds,u 


(ID 


with  the  corresponding  inverse  relation 


Sn(t)  =  \D{,-i;)^dP-dS 

0 


(12) 


=  <7j ,  (0)Z>(0  +  J  D(t  -  &  da^  dt 

o+ 

=  D*dslu 


(13) 


where  the  monotonically  increasing  function  D(t)  is  called  the  uniaxial  creep  compliance  is 
the  measure  of  how  rapidly  the  uniaxial  strain  increases  with  time  under  the  imposed  unit¬ 
sized  step  stress  on=  h(t).  These  relations  are  often  represented  with  the  help  of 
mechanical  analog  models  shown  in  figure  1  as  described  later  on. 
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Figure  1.  Several  analog  models  used  historically  in 
representing  linearly  viscoelastic  material  response. 


1.1.5.  Influence  of  temperature  on  the  viscoelastic  time  scale 

Temperature  is  one  of  the  most  important  environmental  variables  to  affect  polymers  in 
engineering  use,  primarily  because  normal  use  conditions  are  relatively  close  to  the 
material  characteristic  called  the  glass  transition  temperature  -  or  “glass  temperature”  for 
short.  In  parochial  terms  -a  closer  definition  follows-  the  glass  temperature  signifies  the 
thermal  range  in  which  the  material  changes  from  a  stiff  or  hard  material  to  a  soft  or 
compliant  one.  The  major  effect  of  the  temperature,  however  perceived  by  the  user,  is 
through  its  influence  on  the  creep  or  relaxation  time  scale  of  the  material. 

Solids  other  than  polymers  also  possess  characteristic  temperatures,  as,  e.g.  the 
melt  temperature  in  metals;  for  polymers  the  melt  temperature  signifies  melting  of 
crystallites  in  the  (semi-)  crystalline  variants.  Also,  typical  amorphous  solids  like  silicate 
glasses  and  amorphous  metals,  exhibit  distinct  glass  transition  temperatures;  indeed,  much 
of  our  understanding  of  glass  transition  phenomena  in  polymers  originated  in 
understanding  related  phenomena  in  the  context  of  silicate  glasses. 
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1.1. 5.1 


The  entropic  contribution: 


Among  the  long-chain  polymers,  elastomers4  possess  a  molecular  structure  that  comes 
closest  to  our  idealized  understanding  of  molecular  interaction.  The  connection  between 
the  molecular  motion  under  forces  applied  external  to  the  network  can  be  provided  through 
thermodynamic  reasoning.  In  this  way  the  “classical”  constitutive  behavior  under  moderate 
deformations  (up  to  about  100%  strain  in  uniaxial  tension)  has  been  formulated  by  Treloar 
in  19585.  In  the  present  context  the  role  of  the  temperature  is 

of  primary  interest  and  it  suffices,  therefore,  to  quote  his  results  in  the  form  of  the 
constitutive  law  for  an  incompressible  solid  (constant  density).  Of  particular  importance  is 
the  dependence  of  the  stress  on  the  material  property  which  relation  is  for  the  uniaxial 
tension  state  (say,  in  the  1 1 -direction) 


Figure  2:  Deformation  of  a  cube  into  a  parallelepiped.  The  unit  cube  sides  have  been 
stretched  orthogonally  in  length  to  the  extension  or  stretch  ratios  ki,  k2,  and  k3. 


subject  to  A,iA,2A.3=l 


(20) 


where  A,i,  kj  and  A. 3  denote  the  (principal)  extension  ratios  of  the  deformation  illustrated  in 
figure  2  (though  only  for  the  condition  A.i A.2A-3=1),  and  the  multiplicative  factor  consists  of 
the  number  of  chain  segments  between  crosslinks,  N,  Boltzmann’s  constant  k  and  the 
absolute  temperature  T.  Since  for  infinitesimal  deformations  Ai=  1+Sn,  i.e.  the  extension 
ratio  differs  from  unity  by  the  small  axial  strain  sn,  one  finds  upon  expanding  (20)  in 
terms  of  strains  small  relative  to  unity  that  NkT  must  equal  the  elastic  Young’s  modulus 


4  “Elastomer”  is  an  alternative  name  for  rubber,  a  cross-linked  polymer,  possessing  a  glass  transition 
temperature  that  is  considerably  below  normal  environmental  conditions.  Molecule  segments  are  freely 
mobile  relative  to  each  other  except  for  “being  pinned”  at  the  cross-link  sites 

5  Because  this  constitutive  formulation  involves  the  entropy  of  a  deformed  rubber  network,  this  temperature 
effect  of  the  properties  is  usually  called  the  “entropic”  temperature  effect. 
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E,fj .  It  is  thus  established  that  the  (small-strain)  Young’s  modulus  is  directly  proportional 
to  the  absolute  temperature,  and  that  this  holds  also  for  the  shear  modulus  because,  under 
the  restriction/assumption  of  incompressibility  (Xi^2^3=1)  the  shear  modulus  jum  of  the 
rubber  equals  1/3  E,  i.e.  ju^  =  j  E,fj .  Thus  p.fj  / T  =  Nk  is  a  material  constant,  from  which 
it  follows  that  comparative  moduli  obtained  at  temperatures  T  and  T0  are  related  by 

T  T 

/Eol  =—  /EoL  or  equivalently  E,  =  —  E,  (21) 

T  To  l70  lT  To  'T° 

1.1. 5.2  Time-temperature  tradeoff,  the  phenomenon. 

A  much  more  significant  influence  of  temperature  on  the  viscoelastic  behavior  is 
experienced  in  connection  with  the  time  scales  under  relaxation  or  creep.  To  set  the  proper 
stage  we  define  first  the  notion  of  the  glass  transition  temperature  Tg.  To  this  end  consider 
a  measurement  of  the  specific  volume  as  a  function  of  temperature.  Typically,  such 
measurements  are  made  in  a  setting  where  the  temperature  decreases  “slowly”,  since 


Figure  3.  A  typical  set  of  measurements  of  specific  volume  as  a 
function  of  temperature  under  variable  cooling  rates. 


the  rate  of  cooling  has  an  influence  on  the  outcome.  Figure  3  shows  a  typical  result,  which 
illustrates  that  at  sufficiently  low  and  high  temperatures  the  volume  dependence  is  linear, 
with  a  transition  connecting  the  two  segments.  The  glass  transition  temperature  is  defined 
as  the  intersection  of  two  linear  extensions  of  the  two  segments  roughly  in  the  center  of  the 
transition  range.  As  also  indicated  in  that  figure,  an  increase  in  the  rate  of  cooling  causes 
reduced  volume  shrinkage  as  a  result  of  the  unstable  evolution  of  a  molecular 
microstructure  that  consolidates  with  time.  This  phenomenon  is  associated  with  physical 
ageing.  In  practical  terms  the  lowest  -most  stable-  response  curve  is  determined  basically 
by  the  patience  of  the  investigator,  though  substantial  deviations  must  be  measured  in 
terms  of  logarithmic  time  units:  Thus,  relatively  little  may  be  gained  by  reducing  the 
cooling  rate  from  one  to  a  tenth  of  a  degree  centigrade  per  hour. 
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Figure  4:  Idealized  presentation  of  relaxation  behavior  at  various  temperatures  T0  —  T7. 
Scale  markers  represent  a  decade.  The  “experimental  window”  represents  a  typical 
time  range  available  for  a  measurement  sequence. 


We  turn  next  to  the  most  important  effect  of  temperature  on  the  time  scale  and 

present  this  phenomenon  in  terms  of  a  relaxation  response,  say  in  shear.  The  discussion  is 
generic  in  the  sense  that  it  applies,  to  the  best  of  the  collective  scientific  knowledge,  to  all 
time  and  rate  sensitive  properties,  at  least  for  polymers6  For  ease  of  presentation  we  ignore 
first  the  entropic  temperature  effect  discussed  above  and  impose  no  constraint  on  the 
applicable  temperature  range.  Moreover,  we  limit  ourselves  to  considerations  above  the 
glass  transition  temperature,  with  discussion  of  behavior  around  or  below  that  temperature 
range  reserved  for  later  amplification.  This  limitation  is  consistent  with  the  behavior  of  the 
elastomer  considered  in  this  study  as  illustrated  later  on. 

Figure  4  shows  a  series  of  relaxation  curves  measured  at  different  temperatures  in  an  ideal 
test  situation.  Accordingly,  increasing  temperatures  shift  the  relaxation  behavior  to  shorter 
times  parallel  to  the  logarithmic  time  axis.  This  logarithmic  shift  corresponds  thus  to  a 
multiplicative  factor  that  depends  on  the  temperature. 

In  reality  experimental  constraints  do  not  allow  the  full  range  of  relaxation  to  be  measured 
at  any  one  temperature.  Instead,  measurements  can  typically  be  made  only  within  the  time 
frame  of  a  certain  “test  window”  indicated  by  the  dashed  outline  in  figure  4.  Accordingly, 
the  more  realistically  acquired  data  would  appear  as  indicated  in  figure  5.  If  one  accepts 
this  “time -temperature  equivalence”  and  the  associated  “time-temperature  shift”  process, 
then  a  single  curve  — the  master  curve —  may  be  constructed  from  the  segments  in  this 
figure  by  shifting  the  temperature  segments  with  respect  to  one  obtained  at  a  temperature 
that  may  be  chosen  arbitrarily,  so  as  to  construct  the  master  curve  corresponding  to  that 


6  The  technological  evolution  of  metallic  glasses  is  relatively  recent,  so  that  a  limited  amount  of  data  exists 
in  this  regard. 
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Figure  5:  Simulation  of  typical  data  acquired  at  various  temperatures. 


from  which  the  segments  would  have  been  deduced  according  to  the  scheme  in  figure  4. 
Because  this  time-temperature  trade-off  has  been  deduced  from  physical  measurements 
without  the  benefit  of  a  time  scale  of  unlimited  extent,  the  assurance  that  this  “shift- 
process”  is  a  physically  acceptable  or  valid  scheme  can  been  derived  only  from  the  quality 
with  which  the  “shifting”  or  “superposition”  can  be  accomplished.  To  examine  this  quality 
issue  requires  that  test  temperatures  are  chosen  sufficiently  closely,  and  that  the 
measurements  vary  widely  enough  over  the  log-time  range  to  afford  considerable  overlap 
of  the  shifted  curve  segments. 

In  this  construction  a  composite  curve  develops,  tenned  the  “master  curve”,  as  the 
curve  segments  are  shifted  relative  to  one  obtained  at  some  arbitrarily  chosen  temperature. 
At  the  same  time  the  amount  of  shift  along  the  log-time  axis  is  recorded  as  a  function  of 
the  temperature.  This  function,  usually  called  the  “temperature  dependent  shift  factor”,  or 
simply  “shift  factor”  for  short,  is  a  material  characteristic,  and  is  often  designated  by  the 
symbol  .  Figures  7  and  8  illustrate  the  application  of  the  shift  principle  for  a 
polyurethane  elastomer  (Thiokol  Solithane  113,  50/50  composition),  together  with  the 
associated  shift  factor  in  figure  7  with  the  inclusion  of  the  entropic  effect. 


1.1. 5.3  The  role  of  the  entropic  contribution. 

Having  demonstrated  the  “shift-phenomenon”  in  principle,  it  remains  to  address  the 
effect  of  the  entropic  contribution  to  the  time  dependent  master  response.  Recall  that  the 
entropic  considerations  were  derived  in  the  context  of  purely  rubbery  material  behavior, 
and  specifically  in  the  absence  of  viscoelastic  effects.  Thus  any  modulus  variation  with 
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Figure  6:  Relaxation  modulus  for  a  Polyurethane  formulation 
measured  at  various  temperatures  in  uniaxial  tension. 


temperature  is  established,  strictly  speaking,  only  in  the  long  term  time  domain  when 
rubbery  behavior  dominates,  so  that  equations  (21)  apply.  Various  arguments  have  been 
put  forward  to  apply  a  similar  reduction  scheme  to  data  in  the  viscoelastic  transition 
region.  Two  arguments  dominate,  but  they  are  based  on  pragmatic  rather  than  rigorously 
scientific  principles.  The  first  argument  states  that  even  in  the  transition  region  the  polymer 
chain  segments  experience  locally  elastic  behavior  in  accordance  with  the  theory  of  rubber 
elasticity.  Accordingly,  all  curve  segments  obtained  at  the  various  temperatures  should  be 
multiplied  by  their  respective  ratios  of  the  reference  temperature  and  the  test  temperature, 
i.e.  T0/T,  in  the  case  of  modulus  measurements,  and  with  the  inverse  ratio  in  the  case  of 
compliance  measurements,  regardless  of  by  how  many  log-time  units  the  material  behavior 
is  removed  from  the  rubbery  long-tenn  domain.  The  alternative  view  asserts  that  the 
entropic  correction  does  not  apply  in  the  glassy  state  and,  accordingly  should  decrease 
continuously  from  the  long-term,  rubbery  domain  as  the  glassy  or  short-term  behavior  is 
approached.  The  rule  by  which  this  change  occurs  is  not  established  scientifically  either, 
but  is  typically  taken  to  be  linear  with  the  logarithmic  time  scale  throughout  the  transition. 
Ultimately,  the  precision  of  the  data  and  its  treatment  needs  to  decide  as  to  whether  one  or 
the  other  scheme  produces  the  “better  master  curve”  by  the  argument  as  to  which 
procedure  provides  the  best  overlap  of  the  data.  Historically,  the  first  scheme  has  been 
practiced  more  frequently  than  the  latter. 

In  the  present  study  it  turned  out  that  the  material  behavior  was  accessible  mostly 
in  the  near-rubbery  domain,  only  so  that  the  first  of  the  above  reduction  schemes  was 
applied. 
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Figure  7:  Time  temperature  shift  factor  for  reducing  the  Polyurethane  data  in  figure  7  to  that 
in  figure  9.  Tg=  -18°C.  The  solid  line  represents  the  WLF  -  equation 
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Figure  8:  Temperature  reduced  uniaxial  relaxation  modulus  (Master  curve)  for  a 
Polyurethane  formulation  derived  from  data  in  figure  6  and  with 
the  shift  factors  in  figure  7. 

Reference  to  “Range  of  Measurements”  pertains  to  discussion  in  the  text. 


13 


1.1. 5.4 


The  temperature  shift  factor 


While  several  researchers  have  contributed  significantly  to  clarifying  the  concept  and  the 
importance  of  the  time-temperature  superposition  principle,  it  was  the  group  of  Williams, 
Landel  and  Ferry  (1955)  that  has  been  credited  with  formulating  the  time-temperature 
relationship  through  the  now  ubiquitously  quoted  WLF  equation:  they  demonstrated  the 
near  universality  of  this  connection  for  many  diverse  polymers,  and  provided  a  physical 
model  for  the  process  in  terms  of  a  free-volume  interpretation.  Plazek  (1965)  has  supplied 
an  exemplary  demonstration  of  the  nearly  perfect  obeyance  of  the  “shift  phenomenon”  for 
Polystyrene.  Ignoring,  for  brevity  of  presentation,  the  polymer  mechanical  argumentation, 
Williams,  Landel  and  Ferry  (1955)  suggested  that  above  the  glass  transition  temperature 
the  shift  factor  is  given  in  functional  form  by  the  relation 


log[aT] 


Cj  (T  ~  Tfef  ) 

c2  “(T  -Tref  ) 


(22) 


where  Tref  denotes  a  reference  temperature  typically  about  50  °C  above  the  glass  transition 
temperature  of  the  polymer  under  consideration.  The  constants  Ci  and  Ci  vary  from 
polymer  to  polymer,  but  for  many  take  on  values  around  Ci  ~  8.86  and  C2~10 1 .6.  It  is  of 
interest  to  note,  but  of  no  particular  consequence  for  the  purposes  of  this  summary  that  the 
same  functional  relation  can  be  derived  starting  from  thermodynamic  reasoning  without 
any  recourse  to  the  arguments  offered  by  Williams,  Landel  and  Ferry.  As  a  consequence 
one  must  consider,  in  principle,  either  derivation  questionable,  and  the  physical  basis  of 
equation  (22)  as  ultimately  unsettled  and  thus,  as  yet,  empirical. 

An  illustration  of  data  reduction  employing  the  entropic  correction  (applied  over 
the  whole  transition  range)  and  the  time-temperature  superposition  of  measured  data  is 
given  in  figures  6-8.  Figure  6  presents  tensile  relaxation  data  for  a  polyurethane  elastomer, 
possessing  a  glass  transition  temperature  of  -18  °C,  as  measured  at  different  temperatures. 
Shifting  the  curve  segments  relative  to  the  0  °C  data  by  the  shift  factor  shown  in  figure  7 
produces  the  master  curve  in  figure  8.  Because  the  data  segments  were  shifted  relative  to 
the  data  obtained  at  0°C,  the  master  curve  represents  relaxation  behavior  at  that 
temperature.  Data  at  different  temperatures  may  then  be  deduced  with  the  help  of  the  shift 
factor  in  figure  7. 


1.1. 5. 5  Time-Temperature  shifting  near  and  below  the  glass  transition. 

The  time  temperature  shifting  above  the  glass  transition  has  been  presented  as  basically  an 
empirical  rather  than  uniquely  explained  process,  though  many  researchers  firmly  trust  its 
validity  because  of  extensively  consistent  demonstration.  The  applicability  of  the  shift 
principle  to  temperatures  near  and  below  the  glass  transition  has  been  questioned  for  many 
years,  but  is  gradually  gaining  acceptance  with  certain  provisos.  First,  no  functional 
analytic  form  has  been  proposed  -unique  or  conflicting-  that  yields  credence  to  the  effect 
in  terms  of  some  molecular  model. 
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Moreover,  because  phenomena  at  and  below  the  glass  transition  do  not  occur  with 
molecular  conformation  being  in  equilibrium,  the  ideas  underlying  the  shift  phenomenon 
above  the  glass  transition  are  questioned  more  readily,  because  the  role  of  dilatation 
interferes  with  simple  concepts  and  complicates  the  rules  by  which  such  examinations  and 
data  interpretations  are  carried  out.  For  example,  Losi  and  Knauss  have  argued  on  the  basis 
of  free  volume  considerations  that  any  shift  operation  below  the  glass  transition  should 
depend  on  the  temperature  rate  with  which  the  state  of  the  polymer  is  approached,  slower 
rates  leading  to  a  “more  unique”  adherence  to  a  shift  concept.  Many  independent,  though 
not  necessarily  thoroughly  documented  counterparts  have  been  produced  over  the  last 
decade  that  present  equally  supportive  information  of  a  consistent  time-temperature 
superposition  application  through  and  below  the  glass  transition. 


1.1.6  Effect  of  Dilatation  on  Creep  and  Relaxation  Times  (Pressure  effects) 

An  important  consequence  of  free  volume  consideration  is  reflected  in  the  behavior  of 
polymers  under  high  pressures.  In  section  4.4  it  was  noted  that  the  WLF  equation  (22)  can 
be  derived  using  argumentation  based  on  the  role  of  free  volume.  That  reasoning  asserts 
that  an  increase  in  the  space  around  molecule  increases  the  ability  and  speed  of  molecule 
segments  to  adjust  to  new  positions  relative  to  each  other  in  the  process  of  relaxation  or 
creep.  While  this  idea  has  been  considered  mostly  in  connection  with  pressure  superposed 
on  a  tensile  or  shear  deformation,  it  should  be  equally  applicable  to  situations  involving 
positive  dilatation.  This  concept  has  been  investigated  by  only  a  few  investigators,  but 
even  so  an  exhaustive  presentation  of  these  results  is  beyond  the  scope  of  the  present 
objective.  However,  this  concept  plays  a  significant  role  in  the  discussion  of  fracture  and 
stress  wave  propagation.  Stated  as  a  general  thermoelastic  concept  by  Ferry  it  has  been 
explored  experimentally  by  Fillers  and  Tschoegl  (1977),  Moonan  and  Tschoegl  (1985), 
who  studied  the  effect  of  pressure  on  the  relaxation  times  in  the  rubbery  domain,  and  by 
Knauss  and  Emri  (1981,  1987)  and  by  Losi  and  Knauss  (1992)  in  the  formulation  of 
dilatation  induced  nonlinearly  viscoelastic  behavior  for  structural  polymers.  The  most 
recent  confirmation  of  this  phenomenon  or  viscoelastic  feature  is  a  study  by  Qvale  and 
Ravi-Chandar  (2004)  which  addresses  viscoelastic  property  measurements  on  PMMA  and 
PC. 

The  net  result  of  these  investigations  is  that  the  dilation,  whether  positive  or  negative,  has 
an  effect  on  the  material  time  scale  similar  — if  not  identical —  to  a  lowered  or  raised 
temperature,  respectively.  Moreover,  as  is  the  case  of  the  temperature  induced  shift 
phenomenon  for  which  the  free-volume  model  of  the  shift  factor,  equation  (22),  was 
deduced  by  Williams,  Landel  and  Ferry,  the  magnitude  of  the  dilatation  necessary  to  effect 
significant  changes  in  the  relaxation  or  creep  time  scales  may  be  very  small,  i.e.,  on  the 
order  of  only  10'3  or  less.  With  regard  to  failure  and  /or  shock  wave  environments  this  is  an 
important  observation  because  at  the  tip  of  a  crack  the  high  stresses  are  thus  automatically 
a  source  of  shift  in  time  response  in  the  material.  This  notion  appears  to  be  particularly 
important  for  the  fracture  of  structural  polymers  because  the  high  crack  tip  stresses 
dominate  and  accelerate  the  viscoelastic  yield-like  behavior  of  the  materials  in  the  region 
where  fracture  occurs.  With  respect  to  the  propagation  of  pressure  waves  this  notion 


15 


implies  an  automatic  shift  of  the  relaxation  or  creep  spectrum  to  the  more  dissipative 
region  of  the  spectrum  if  an  elastomer  is  involved  such  as  is  the  topic  of  this  investigation. 
Associated  with  this  shift  to  higher  dissipation  an  apparent  stiffening  of  the  material  is 
experienced. 


1.1. 7  Representation  of  the  Relaxation  and  Creep  Functions  (Spectra  and  functional 
representation) 

Various  mathematical  forms  have  been  suggested  and  used  to  represent  the  material 
property  functions.  Preferred  forms  have  evolved  with  precision  being  balanced  against 
ease  of  mathematical  use  or  minimum  number  of  parameters  required.  All  viscoelastic 
material  functions  possess  the  common  characteristic  that  they  vary  monotonically  with 
time:  relaxation  functions  decrease  and  creep  functions  increase  monotonically.  A  second 
characteristic  of  realistic  material  behavior  is  that  time  is  (almost)  invariably  measured  in 
terms  of  (base  10)  logarithmic  units  of  time. 

Early  representations  of  viscoelastic  responses  were  closely  allied  with  (simple) 
mechanical  analog  models  (Kelvin,  Voigt)  or  their  derivatives.  Without  delving  into  the 
details  of  this  evolutionary  process,  their  generalization  to  broader  time  frames  lead  to  the 
“spectral  representation”  of  viscoelastic  properties,  so  that  it  is  useful  to  present  only  the 
rudiments  of  that  development.  The  “building  blocks”  of  the  analog  models7  are  the 
Maxwell  and  the  Voigt  models  illustrated  in  figure  la  and  lb.  Under  a  step-wise  applied 
deformation  of  magnitude  s0  — represented  by  the  separation  of  the  force-application 
points  in  the  Maxwell  model,  Figure  la —  the  stress  (force)  abates  or  relaxes  by  the 
relation 


r(t)  =  £0  •  exp(- 


(14) 


where  |  =  ^m/  is  called  the  relaxation  time.  Similarly,  applying  a  step  stress 

/  7  m 

(force)  of  magnitude  t0  to  the  Voigt  element  (Figure  lb)  engenders  a  time  dependent 
separation  (strain)  of  the  force-application  points  described  by 


£(0  =  —  [1-exp  -  Vj] 


(15) 


where  £  = 


is  now  called  the  retardation  time  since  it  governs  the  rate  of  retarded  or 


delayed  motion.  Note  that  this  representation  is  symbolic  for  illustration  purposes  here  and 
that  the  retardation  time  for  the  Voigt  material  is  not  necessarily  meant  to  be  equal  to  the 
relaxation  time  of  the  Maxwell  solid.  In  fact,  in  general  relaxation  times  are  not  equal  to 


7  Force  F  ~  stress  x  and  deformation  5  -strain  s. 
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retardation  times.  The  relaxation  modulus  and  creep  compliance  commensurate  with 
equations  14  (Maxwell  model)  and  15  (Voigt  model)  for  the  Wiechert  and  Kelvin  models 
(Figures  lc  and  Id)  are,  respectively 


Kt)  =  Aoo  +  Z  Mn  eXP(_/£;) 


and 


J(t)  —  J  a  +Z- 


Mn 


l-exp(-  /  ) 

/  ?  n 


+  rj0t. 


(16) 

(17) 


Here  Jg  and  q„  arise  from  letting  rj  i — ; >0  (the  first  Voigt  element  degenerates  to  a  spring) 
and  pn— >0  (the  last  Voigt  element  degenerates  to  a  dashpot).  These  series  representations 
with  exponentials  are  often  referred  to  as  Prony  series. 

As  the  number  of  relaxation  times  increases  indefinitely,  the  generalization  of  the 
expression  for  the  shear  relaxation  modulus,  becomes 

00  /  d<5 

A(t)=  Jtf(£)exp(-  jA-f  (18) 

o  £ 


where  the  function  H(g)  is  called  the  “distribution  function  of  the  relaxation  times”,  or 
relaxation  spectrum,  for  short;  the  creep  counterpart  presents  itself  with  the  help  of  the 
retardation  spectrum  £(Q  as 


J(t)  =  Jg  +  J  L  (q)[  1  -  exp(-  Zg) ]  %  +  1!  ■  t . 


(19) 


d£ 

o  -  -  C 

Note  that  although  the  relaxation  times  £,  and  the  retardation  times  C  do  not,  strictly 
speaking,  extend  over  the  range  from  zero  to  infinity,  the  integration  limits  are  so  assigned 
for  convenience  since  the  functions  H(q)  and  L(Q  can  always  be  chosen  to  be  zero  in  the 
corresponding  part  of  the  infinite  range. 


A  typical  example  of  a  relaxation  modulus,  say  in  shear  (or  in  uniaxial  tension)  is 
illustrated  in  figure  8.  It  is  particularly  noteworthy  that  the  time  scale  ranges  over  many 
decades  which  are  typically  very  difficult  or  impossible  to  access  directly  in  a  standard 
laboratory  environment.  It  is  the  subject  of  this  report  to  display  this  complete  time  scale 
on  the  basis  of  the  Time-Temperature  Tradeoff  principle,  as  discussed  below.  While  we 
have  chosen  to  perfonn  this  characterization  in  tenns  of  the  relaxation  phenomenon,  it 
should  be  stated  that  other  laboratories  might  perform  this  task  in  tenns  of  so-called 
dynamic  or  complex  environment.  This  mode  is  not  discussed  here,  though  interrelations 
can  be  derived  immediately  from  equations  (1-13). 
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Figure  9:  Example  of  a  typical  relaxation  modulus  function. 


1.1.8.  General  Three-Dimensional  Constitutive  Description. 

In  sections  2. 7.2. 3  through  2. 7. 2. 5  we  have  dealt  with  one-dimensional  relations  between 
stress  and  strain  components.  We  now  generalize  these  descriptions  for  fully  three- 
dimensional  interactions  between  stresses  and  strains.  We  do  this  here  only  for  isotropic 
material  description  and  observe  that  in  the  limit  of  vanishing  viscoelasticity,  i.e.  for  purely 
elastic  behavior,  the  constitutive  laws  need  to  revert  to  their  elastic  counterparts.  This 
representation,  as  illustrated  so  far  for  one-dimensional  descriptions,  essentially  replaces 
the  elastic  constants  by  the  convolution  operators.  One  notes  further  that  elastic  isotropy 
requires  two  independent  material  constants  for  a  complete  material  characterization. 
Because  in  the  limit  of  elastic  behavior  the  viscoelastic  response  must  revert  to  elastic 
behavior,  it  stands  to  reason  that  isotropic  viscoelasticity  requires  two  independent  material 

o 

functions  for  a  full  characterization  .  Formulating  the  full  constitutive  response  is  most 
readily  accomplished  in  terms  of  the  deviatoric  (shear)  and  bulk  (volumetric)  property 
functions  pit)  and  k(t).  Bearing  in  mind  that  for  three-dimensional  problems  the  stress  and 
strain  are  function  of  both  space  x  (=xi,X2,X3)  and  time  t,  the  isotropic  linearly  viscoelastic 
solid  constitutive  formulation  is  then,  in  terms  of  the  deviatoric  and  dilatational 
representation, 

r  Se(x,c) 

Sy(x,t)  =  f u  h  dg  (23) 

0 


This  argument  clearly  does  not  constitute  a  proof. 
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(24) 


r„(M)  =  3ft(<'-f)a[e(X^)~aA7'(X-^ 

0  ^ 

where  use  has  been  made  of  the  general,  temperature  reduced  time  t’  and  a  time-  or 
rate-independent  volumetric  coefficient  of  thermal  expansion  a  has  been  assumed  for  ease 
of  presentation.  If  a  history-dependent  coefficient  of  thermal  expansion  is  involved  then 
aAT  must  be  replaced  by  a  convolution  equation.  Because  the  integration  variable  is  2  and 
not  4’  which  depends  on  these  integrals  are  not  Duhamel  or  convolution  integrals. 
Nevertheless,  we  shall  use  the  short  notation  of  the  star  operation  and  write,  in  the  absence 
of  thermal  expansions, 


Sij(x,t)  =  ju*deij(x,£)  (25) 

and  Ta(x,f)  =  3 k*  dskk  (x,  g).  (26) 

Alternate  forms  in  terms  of  the  non-deviatoric  stress  and  strain  components  are  then, 
without  the  thermal  expansion  components 

T](x,t)  =  2ju*dci/x,  g)  +  dyi k— 2/3  \ifdckkfx,  g)  (27) 

where  fi(t)  and  k(t)  represent  the  shear  and  the  bulk  moduli  described  in  sections  2. 1  — 2. 

Equivalently,  one  can  define  the  constitutive  behavior  in  terms  of  the  Lame  form  if  the 
elastic  Lame  parameter  2  is  considered  to  be  a  viscoelastic  function  so  that 


Tij(x,t)  =  X*dskk(x,g)  dy  +  2 \i*dsij{x,t).  (28) 

Interconnections  between  A,(t)  and  the  other  moduli  may  be  deduced,  but  are  not  of 
immediate  consequence. 


2.  MEASUREMENTS  OF  ELASTOMER  PROPERTIES 

Having  reviewed  the  general  principles  of  linear  viscoelasticity,  we  turn  next  to  the 
determination  of  the  viscoelastic  characteristics  of  the  elastomer  at  hand.  Following  an 
analysis  of  the  measurement  precision  for  reasons  of  demonstrating  the  reliability  of  the 
measured  functions,  we  document  the  master  curve(s)  for  the  elastomer,  along  with  a 
detailed  description  fo  the  thermo-mechanical  trade-off  between  time  and  temperature,  so 
that  application  of  the  data  can  be  made  in  any  thermal  environment. 


2. 1  Tools. 

The  viscoelastic  characterization  was  accomplished  through  measuring  the  relaxation 
behavior  in  uniaxial  extension.  Such  measurements  require  the  recording  of  the  force  on  a 


19 


tensile  coupon  as  a  function  of  time  resulting  from  a  “rapidly”  applied  strain  that  is  held 
constant  with  time.  Such  measurements  are  in  principle  readily  accomplished  with  the  aid 
of  a  (servo-)  hydraulic  tensile  tester  (MTS).  Though  alternate  straining  devices  can  be 
constructed  for  greater  precision  in  strain  definition,  this  mode  was  chosen  to  minimize  the 
time  and  expense  for  the  overall  test  program. 


2.1.1  Calibration  of  load  cell 

The  tests  were  performed  on  an  MTS  load  unit  Model  #  358.1  with  an  axial  load 
transducer  by  Interface,  Model  #  1500  ASK-200  having  a  capacity  of  200  lbf.  The 
controller  was  the  installed  TestStar  Tim  system  relying  on  Model  793.00  system  software. 

Prior  to  starting  test  the  load  cell  was  calibrated  statically  with  calibrated  weights.  Since 
the  load  range  of  the  load  cell  extended  to  200  lbs  and  the  lowest  load  measured  in  the 
relaxation  tests  was  on  the  order  of  15  lbs,  tensile  loads  could  be  determined  to  a  fraction 
of  a  percent,  which  was  consistent  with  the  error  adhering  to  the  determination  of  the 
strains  (see  later). 


Caliper  Calibration 


MTS  indication,  mm 


Figure  10:  Calibration  of  the  crosshead  motion  of  the  MTS  machine 
Against  the  gage  strain  measured  with  calipers. 


2. 1. 2  Test  Strain  Magnitude 

Because  it  was  anticipated  initially  that  measurements  would  be  conducted  into  the  glassy 
domain  for  which  ancillary  studies  by  Lu  and  Knauss  (1999),  and  by  Knauss  and  Zhu 
(2002,  2003)  had  shown  that  nonlinear  behavior  can  be  encountered  at  strains  of  only  half 
a  percent,  it  was  desired  to  employ  strains  on  that  order.  However,  it  turned  out  that  the 
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precision  of  the  servo-controller  of  the  test  machine  was  smaller  than  stated  and 
anticipated.  Many  tests  had  been  perfonned  with  spurious  results  which  at  times  indicated 
that  physical  properties  changed  with  time  and  under  repeat  use  in  an  apparently 
systematic  manner.  It  was  only  after  a  careful  calibration  of  this  control  capability  and  after 
numerous  uncertain  and  spurious  variations  that  the  error  was  cleared  up.  As  a 
consequence  of  this  finding,  the  strain  amplitude  for  the  relaxation  tests  needed  to  be 
increased.  However,  it  was  judged  that  the  penalty  of  exceeding  the  possibly  nonlinear 
range  was  minor  compared  to  the  need  for  sufficient  precision.  The  penalty  concern  was 
largely  eliminated  when  it  was  found  at  the  end  of  the  measurements  that  the  elastomer 
never  entered  into  the  glassy  domain,  but  remained  in  the  near  rubbery  region  where  the 
range  of  linear  behavior  is  significantly  larger  than  in  the  glassy  domain.  After  the  change 
to  larger  strain  values  was  made,  all  measurements  became  systematic  and  in  accordance 
with  expected  viscoelastic  behavior. 


2.1.3  Precision  of  strain  prescription  (strain  vs.  displacement  monitoring) 

It  was  desired  to  determine  both  the  axial  and  the  lateral  strain  in  the  specimens  under  load. 
To  this  end  the  Digital  Image  Correlation  (DIC)  method  was  employed,  specifically  so  as 
to  be  able  to  estimate  the  lateral  strain  (change).  It  ultimately  turned  out  that  the  precision 
of  that  method  at  the  strain  amplitudes  used  was  insufficient  to  allow  an  estimation  of  a 
more  than  nominal  Poisson  effect  to  be  determined.  Thus  while  the  data  of  this  quantity 
was  recorded  as  matter  of  routine,  the  result  is  essentially  that  a  (near-)  incompressible 
behavior  was  found  (Poisson’s  ratio  v~0.5).  For  the  axial  strain  determination  the  use  of 
the  servo-hydraulic  test  machine  required  really  prescription  of  the  length  change  of  the 
specimen.  To  this  end  a  test  series  was  performed  in  which  the  routinely  and  consistently 
gripped  specimen  would  be  extended  by  recorded  amounts.  The  strain  in  the  specimen  test 
section  was  recorded  by  measuring  the  separation  of  scribe  marks.  In  addition,  the  axial 
strain  was  determined  by  the  DIC  method,  and  the  strain  was  plotted  as  a  function  of  the 
test  ram  displacement  as  indicated  on  the  instruments  monitoring  the  test  machine.  Figure 
10  shows  the  relation  between  the  test  frame  supplied  displacement  and  the  test  section 
strain.  This  relation  was  used  as  a  guide  in  configuring  tests,  yet,  ultimately,  the  directly 
measured  strain  (via  DIC)  was  used  for  data  reduction. 

A  note  is  in  order  on  the  method  employed  here  for  the  determination  of  the  strains  by 
means  of  DIC.  In  principle  DIC  is  able  to  determine  simultaneously  displacements  as  well 
as  (two-dimensional)  strain  fields.  This  capability  of  DIC  is  welcome  when  inhomogenous 
strain  fields  prevail.  However,  when  the  strain  field  is  homogeneous  it  has  turned  out  in 
our  past  experience  that  the  precision  of  the  strain  is  increased  if  the  DIC  deduced 
displacements  are  processed  separately  for  calculating  the  strain,  at  the  expense  of 
considerable  additional  time  investment.  However,  it  was  felt  that  for  this  program  it  was 
desirable  to  maximize  precision  over  quantity  or  speed  of  data  acquisition,  so  that  the  more 
intensive  and  slower  process  was  chosen. 
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Figure  11:  Roster  of  the  centers  of  image  regions  within  which 
displacements  were  determined  by  means  of  the  DIC  method. 


To  describe  the  strain  determination  refer  to  Figure  1 1  in  which  numbers  along  the 
ordinate  and  along  the  abscissa  represent  the  centers  of  small  domains  or  image  regions 
around  which  displacements  (and  strains)  were  determined  by  the  DIC  method.  Because 
the  deformations  in  a  uniaxial  strain  situations  are  homogeneous  the  y-displacements  (v) 
were  averaged  across  the  roster  (left  to  right:  m=T-toM)  for  each  line  n=l-to-N  parallel  to 
the  abscissa;  subsequently,  these  x-averaged  values  were  plotted  (Excel)  against  the 
ordinate  number  and  this  essentially  linear  relation  was  then  fitted  with  a  least-square 
straight  line  (also  Excel)  with  the  slope  being  the  desired  strain.  An  example  plot  of  this 
type  is  shown  in  Figure  12.  We  have  also  explored  averaging  the  DIC  determined  strains 
instead  of  proceeding  via  the  displacements,  yet  determined  that  the  chosen  method  was 
somewhat  advantageous  over  the  averaging  method. 


2.1.4  DIC  image  file  size  and  strain  resolution 

In  an  attempt  to  improve  the  resolution  of  the  DIC  method  for  an  improved  acquisition  of 
the  lateral  strains  a  camera  of  higher  pixel  count  was  used  for  comparison  purposes.  To 
this  end  several  comparison  measurements  were  executed  with  a  3.3  Megabite  camera 
(Nikon  Coolpix  995)  and  a  5.4  Megabite  counterpart  (Nikon  DX1).  One  might  expect  from 
first  principles  that  a  higher  pixel-per-mm  count  should  lead  to  a  higher  precision  in  the 
strain  resolution,  yet  such  was  not  clearly  the  case  by  the  DIC  method  through  which 
strains  were  determined.  The  demonstration  of  that  fact  is  illustrated  in  Figure  13.  Whether 
the  strain  accuracy  were  to  improve  through  the  use  of  higher  pixel  density  in  the  range 
available  if  the  DIC  method  was  used  directly  to  deduce  the  strain  was  not  explored,  since 
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the  present  “optimal”  method  yielded  no  improvement  for  the  strain  extraction. 


485-487-5MM-0C 


360  380  400  420  440  460  480 


y-Position,  pixels 

Figure  12:  Illustration  of  the  data  fit  from  the  DIC  method 
for  determining  the  axial  strain  (  =3.25  %  in  this  example). 


For  the  record  is  should  be  stated  that  the  strains  for  the  two  pixel  density  records  was 
subjected  to  further  analysis  by  fitting  least-square  straight  lines  through  the  data 
incorporating  the  origin.  The  3.3  Mpixel  data  was  represented  by  a  proportionality  between 
grip  displacement  and  strain  with  a  factor  of  0.649  (R2=0.996)  and  for  the  5.2  Mpixel 
count  by  0.669  (R  =0.988);  these  fits  represented  a  standard  deviation  from  the  straight 
lines  of  0.99  for  the  lower  pixel  density  and  0. 188  for  the  higher  one.  Thus,  the  precision  of 
the  strains  determined  either  way  are  precise  to  within  one  or  two  tenths  of  a  percent. 
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Figure  13:  Comparison  for  DIC  deduced  strains  with  the  aid 
of  different  camera  image  resolutions. 


2.1.5  Temperature  control  and  tracking 

Temperature  control  was  achieved  via  a  Russel’s  environmental  chamber  (Model  #  RDB- 
3-LN2-.33)  employing  a  MICRISTAR  controller  (Model  #828-D10-403-901-130-00). 

To  reach  temperatures  below  0°  C  liquefied  nitrogen  was  needed.  At  temperatures  below  - 
55  to  60  °C  this  cooling  method  became  troublesome  and  rather  expensive  because  of  the 
large  amount  of  liquid  Nitrogen  used. 

A  TestStar  controller  was  used  in  conjunction  with  an  OMEGA  CN76000  Microprocessor- 
based  Temperature/process  controller  and  a  type  J  Thermocouple  located  near  the 
specimen  inside  the  chamber.  At  low  temperatures  window  frosting  was  a  standard 
occurrence,  but  a  fan  heater  kept  the  multi-pane  window  clear  for  optical  (DIC-camera) 
access. 

It  was  typically  difficult  with  the  existing  equipment  to  predetermine  the  temperature 
precisely.  This  fact  resulted  from  the  restriction  of  the  temperature  controller  and  the  long 
time  scale  required  to  achieve  a  prescribed  temperature.  Tests  were  thus  conducted  to 
assure  a  stable  test  temperature.  These  test  runs  established  that  a  rather  long  time  (hours) 
were  required  to  have  the  temperature  control  chamber  reach  a  constant  temperature  near 
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the  desired  value,  and  that  the  difference  between  the  pre-set  (“dialed-in”)  value  and  the 
finally  achieved  one  was  not  the  same.  Inasmuch  as  the  precise  pre-detennination  of  the 
temperature  was  more  a  matter  of  convenience  than  a  scientific  constraint,  it  was  decided 
that  it  was  more  important  to  reach  a  steady  temperature  rather  than  a  particular, 
predetennined  value.  Correspondingly,  the  temperature  was  monitored  continually  in  two 
ways:  One  was  by  means  of  the  temperature  record  built  into  the  test  chamber,  with  the 
thennocouple  near  the  entrance  of  the  cooled  or  heated  air.  The  other  derived  from  an 
additional  thennocouple  placed  close  to  the  test  specimen,  the  output  of  which  was 
monitored  continuously  as  a  function  of  time.  The  two  methods  resulted  in  close 
agreement,  and  thus  the  output  of  the  thennocouple  that  was  placed  close  to  the  specimen 
was  used  for  thennal  definitions.  As  a  result  of  this  process  the  test  results  are  reported  for 
temperatures  that  are  rarely  integral  values  of  temperature,  but  typically  include  fractional 
values 

Temperatures  were  tracked  continuously.  In  the  early  stages  of  the  investigation  relaxation 
data  was  rejected  if  the  temperature  variation  was  judged  too  large  -in  excess  of  +/-0.5°C. 
As  test  progressed  that  situation  occurred  less  and  less,  so  that  the  final  data  is  essentially 
free  of  concerns  in  that  direction.  On  the  whole,  analysis  of  the  temperature  history 
established  that  the  standard  deviation  for  temperatures  in  the  steady  state  (long  term  tests) 
were  on  the  order  of  1/3  °C. 


2.2  Properties  Determination 

To  minimize  untoward  and  poorly  understood  nonlinear  response  an  effort  was  made  to 
prescribe  sufficiently  small  deformations  for  the  test  sequences.  This  requirement 
supported  the  choice  of  large  specimens  so  as  to  be  able  to  prescribe  strains  accurately 
through  machine  controlled  grip  displacements.  However,  because  failure  would  occur 
typically  under  “large”  strains,  a  different  specimen  configuration  was  required  that  started 
with  a  shorter  test  section,  so  that  the  stroke  of  the  machine  was  sufficient  to  impose  failure 
strains  at  all  temperatures  considered. 


cross  sectional  dimensions  each  specimen  was  measured  separately. 
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2.2.1  Specimen  configurations  for  relaxation  and  rupture  tests. 

The  test  specimen  is  illustrated  in  figure  14.  Specimens  were  machined  from  nominally 
two  mm  thick  sheet  stock  on  an  NC-milling  machine  with  a  0.5  inch  diameter,  helically 
fluted  cutter  at  3600  rpm.  The  specimen  was  held  in  a  fixture  which  guaranteed  proper 
backing  against  the  cutting  force.  The  dimensions  of  each  specimen  were  measured 
individually  for  properties  evaluations.  Specimens  were  used  for  measuring  relaxation  at 
different  temperatures,  but  not  more  often  than  for  five  tests.  In  an  initial  evaluation  aimed 
at  ascertaining  whether  residual  deformations  or  property  changes  were  encountered  upon 
repeated  use  gave  erroneous  messages,  because  the  test  machine  provided  inaccurate 
control  (see  elsewhere  in  this  report).  Upon  correcting  for  this  deficiency  it  was  established 
that  multiple  use  of  the  specimens  was  permissible  at  the  strain  levels  encountered. 


30  typ 


Figure  15:  Specimen  dimensions  for  failure  tests,  in  mm,+/-  0.3  mm. 
Each  specimen  was  measured  separately  for  each  test. 


2.2.2  Loading  (straining)  programs  for  relaxation 

The  MTS  machine  was  programmed  to  strain  the  specimens  in  a  ramp  fashion  ,  with  the 
constant  strain  (displacement)  being  achieved  in  about  a  second.  As  a  consequence  of  this 
time  scale  the  data  between  1  and  10  seconds  was  not  considered,  since  it  was  too  heavily 
flavored  by  the  loading  transient.  Moreover,  it  was  not  deemed  necessary  to  employ 
special  methods  to  extract  relaxation  data  from  this  short-time  history,  since  the  tests  could 
be  conducted  for  about  a  day  each,  thus  rendering  data  in  excess  of  three  decades  without 
even  the  first  decade  influenced  by  the  loading  transient. 
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2.2.4  Relaxation  at  various  temperatures 


Relaxation  tests  were  thus  conducted  at  various  temperatures.  The  results  are  summarized 
in  figure  16.  Note  that  the  ordinate  represents  the  temperature -reduced  scale,  in  preparation 
for  the  next  step  in  the  data  reduction. 


2.2. 5  Master  relaxation  curve  and  time-temperature  shift  data 


The  data  from  figure  17  has  been  shifted  according  to  the  time -temperature  shift  principle. 
The  result  is  shown  in  figure  18,  along  with  the  resulting  shift  factor  in  figures  19  and  20. 
We  note  first  that  the  superposition  fit  of  the  temperature  segments  in  figure  17  render  a 
rather  smoothly  meshing  set  of  curves  that  produces  a  distinct  master  curve.  It  is  this 
quality  of  fitting/meshing  that  typically  assures  that  applicability  of  the  shift  process.  Next, 
note  that  the  shift  factor  data  in  figure  19  is  represented  well  by  the  (standard)  WLF 

equation  given  by  (0  =  temperature,0  C ) 


log  aT  = 


8.86(0-0.6) 
101.6 +  (0-0.6)  ‘ 


(29) 


Williams,  Landel  and  Ferry  established  for  about  50  polymers  that  the  glass  transition 
temperature  is  50°C  below  the  reference  temperature,  which  is  established  here  as  -0.6  °C. 
I  ti  sclear  by  consulting  figure  xxx,  that  a  variation  of  +/-  5  °C  in  the  glass  transition  yields 
a  readily  distinguishable  difference  in  the  data  fit.  This  value  may  differ  by  a  small  amount 
from  a  determination  by  differential  scanning  calorimetry  (DSC),  though  that  method  has 
its  own  problems  with  being  more  precise. 


2.2.6  Poisson  Behavior 

To  deduce  an  estimate  of  the  dilatational  behavior  an  attempt  was  made  to  measure  the 
time-dependent  Poisson  effect.  It  turned  out  that  the  DIC  method  was  insufficiently 
accurate  to  a  definitive  statement  about  this  behavior.  While  some  test  sequences  resulted 
in  the  expected  behavior  of  increasing  Poisson  number  with  time  under  relaxation,  there 
were  other  tests  in  which  the  opposite  or  a  constant  behavior  was  observed.  All 
measurements  were  close  to  (nearly)  incompressible  behavior  with  a  contraction  ratio  of 
close  to  14  being  observed.  This  kind  of  measurements  is  notoriously  difficult,  especially  in 
the  long-time  or  rubbery  domain,  where  very  small  deviations  from  a  value  of  14  leads  to 
large  variations  in  the  bulk  response. 
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Figure  16:  Relaxation  response  at  various  temperatures 
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Figure  17:  Master  Relaxation  Curve,  referred  to  0°C  (0.22°C) 
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Figure  18.  Shift  factor  data  and  fitted  WLF  equation  (see  text) 
The  glass  transition  temperature  is  about  -50.6  °C 


Temperature  C 


Figure  19:  Shift  Factor  data  compared  with  the  WLF  equation  for 
three  values  of  the  glass  transition  temperature 


2.2.7 


Rupture  Behavior 


An  important  characteristic  in  engineering  failure  analyses  are  “properties”  determined  to 
assess  the  stress  or  strain  capability  of  a  material.  The  purpose  of  such  measurements  is 
typically  to  bound  the  capability  of  the  material  or  to  establish  “allowable”  quantities. 
Most  often  this  is  effected  through  uniaxial  tensile  testing,  though  compression  testing  is 
also  of  importance  in  appropriate  engineering  situations.  Just  to  illustrate  the  effect  of  rate 
of  deformation  (and  temperature  by  way  of  the  time-temperature-tradeoff  Principle)  on 
such  allowables  as  the  ultimate  stress  and/or  strain  achieved  in  uniaxiale  tensile  tests 
under  constant  rate  of  deformation  are  presented  here.  Considerably  different  values  are 
achieved  depending  on  the  time  scale  involved  in  the  measurements  (rate  ~  f '). 


2.2. 6. 1  Definition  of  an  equivalent  constant  strain  rate. 

Failure  strain  determination  requires  a  specimen  (cf.  figure  15)  that  allows  failure  to 
occur  away  from  stress  concentration  caused  by  the  grips  of  the  test  machine.  The 
resulting  inhomogeneity  of  the  specimen  straining  leads  to  a  (large)  strain  that  is  not 
proportional  to  the  grip  motion.  Figure  20  shows  a  typical  example  of  a  test  result,  with 
the  measured  values  of  the  strain  across  benchmarks  denoted  on  the  specimens.  These 
measurements  were  conducted  near  room  temperature  (23  °C)  where  relaxation  plays  a 
minor  role,  though  the  latter  has  no  influence  on  the  geometric  effects  incorporated  the 
results  shown  in  figure  21.  Clearly  the  behavior  is  anything  but  linear.  To  establish  a 
quasi-constant  strain  rate,  we  replace  the  true  strain  curve  by  a  least  square  fit  by  a 
constant-slope  line  to  define  a  (constant)  strain  rate.  This  is  roughly  justified  inasmuch  as 
the  viscoelastic  response  is  really  sensitive  to  the  logarithm  of  the  strain  rates,  rather 
than  relatively  small  variations  of  that  quantity. 

The  next  figure  shows  a  typical  trace  of  the  nominal  (engineering)  stress  as  a  function  of 
the  applied  strain.  The  rise  in  the  stress  at  elevated  strain  (-200  to  400  %)  arises  from  the 
molecular  orientation  typical  of  high  strain  elastomer  behavior. 


2. 2. 6. 2  Results  for  uniaxial  failure. 

Two  measures  for  failure  were  considered,  namely  ultimate  strain  and  ultimate  stress.  In 
the  following  figures  strain  is  given  in  terms  of  engineering  values,  i.e.  (l-ln)/l0.  and  not  in 
percent  values  Stresses  are  given  in  terms  of  the  engineering  stress  as  well  as  in  terms  of 
true  stress,  which  was  deduced  from  the  extension  and  the  lateral  contraction  based  on 
volume  constancy.  This  was  used  regardless  of  the  test  temperature.  The  strain  rate 
(abscissa)  is  computed  for  the  three  nominal  strain  rates  at  different  temperatures. 
However,  the  strain  rates  have  been  converted  to  effective  values  according  to  the  time- 
temperature  trade-off  principle  represented  in  figure  18. 
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One  notes  that  the  maximum  strain  (Figure  24)  remains  on  the  order  of  3  (300%),  though 
there  is  a  small  trend  to  decreased  strain  capability  as  higher  effective  strain  rates  are 
encountered. 

The  stresses,  on  the  other  hand,  appear  to  pass  through  a  maximum  at  a  strain  rate  of  104 
per  sec,  with  a  local  minimum  at  106  per  sec.  this  holds  true  for  the  engineering  and  for 
the  true  stress.  The  three  strain  rates  used  in  the  tests  give  rise  to  a  systematic  shift  to 
higher  rates  as  one  passes  from  5  to  50  mm/sec  rates.  However,  the  stress  amplitudes 
vary  considerably;  experience  would  indicate  that  this  result  derives,  in  good  portion, 
from  micro  fractures  incurred  in  the  specimen  machining  process. 

As  a  minimum,  these  data  provide  a  set  for  response  bounds  on  the  material  in 
engineering  applications. 
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Figure  20:  Definition  of  an  equivalent  strain  rate. 
Strain  measured  in  the  test  section  of  the  specimen. 


2.2. 8  The  influence  of  Pressure  on  time  dependent  behavior 

In  section  1.1.6  it  was  observed  that  pressure  has  a  notable  influence  on  the  relaxation 
times  of  a  viscoelastic  polymer,  similar  to  the  influence  of  temperature:  Higher  pressure 
changes  the  internal  time  scale  (relaxation  times)  to  longer  times  similar  to  cooling  the 
material.  The  only  direct  experimental  study  of  this  effect  has  been  supplied  by  Fillers 
and  Tschoegl  (1977),  and  by  Moonan  and  Tschoegl  (1985). 
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Figure  21:  Example  of  an  engineering  (nominal)  stress-strain  response.  Note  the 
stiffening  of  the  material  at  higher  strains  (molecular  orientation  and  toughening) 


The  effect  of  changing  pressure  and  temperature  is  included  through  the  concept  to  the 
intrinsic  material  time  defined,  for  a  change  in  pressure  and  temperature  to  other  constant 
values  by 


t 


t'= 


(30) 


a 


If  changes  in  pressure  and/or  temperature  occur  continuously,  this  relation  applies 
instantaneously  so  that  as  a  function  of  any  pressure  or  temperature  history  one  finds  the 
pressure/temperature  reduced  time  from 


(31) 


0  ttlp  p 

In  these  equations  the  shift  factor  is  given  by 


a 


T,P 


aTP[T{t),P{t)\ 


(32) 


and  in  detail  by 


(33) 


where  the  function  0(P)  is 
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Figure  22:  Rupture  engineering  stress  as  a  function  of  the  rate  of  deformation 
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log  10  (strain  rate  x  shift  factor) 

Figure  23:  True  rupture  stress  as  a  function  of  the  rate  of  deformation. 
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Figure  24:  Rupture  strain  as  a  function  of  deformation  rate 


(34) 


_,n,  .  _  1  +  cAP  _  .  r  1  +  CcP  _ 

0{P)  =  c3  ln[- - 4-]  -  c5  ln[—  4_] 


1  +  c4P0 


1  +  C6^0 


where  P  is  measured  in  kbars  and  the  constants  c<  are,  for  example,  for  Hypalon  40,  one 
of  the  materials  studied  by  Fillers  and  Tschoegl  (1977),  Ci=2.66,  C2=60.0,  03=152.0, 
04=0.526,  c5=105.8,  c6=0.576;  similarly  it  was  found  that  ci=2.24,  00=59.4,  03=145.6, 
C4=0.402,  05=69.4,  C6=0.395  for  Viton  B.  For  more  detailed  information  the  reader  is 
referred  to  the  publicaitons  by  Tschoegl  and  his  coworkers. 


It  bears  mentioning  that  since  C4  and  Ce  are  less  than  unity  so  that  for  small  pressures  the 
lat  relation  is  basically  linear  in  the  pressure.  However  once  pressures  exceed  several 
kbars  that  is  no  longer  the  case,  though  the  nonlinearity  is  “mild”. 


2.2. 8. 1  Examples  of  pressure  shifting 

Consider  three  examples  to  show  the  effect  of  pressure  in  tenns  of  the  isothermal 
pressure  shift  factor  in  figures  25,26  and  27. 


Figure  25:  Pressure  shift  factor  for  Hypalon  40. 
(Reference  temperature  25  °C,  and  pressure  at  1  bar 

Fillers  and  Tschoegl,  1977) 
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Figure  26:  Pressure  shift  factor  for  Viton  B. 
(Reference  temperature  25  °C,  and  pressure  at  1  bar 

Fillers  and  Tschoegl,  1977) 


Figure  27:  Pressure  shift  factor  for  Neoprene  WB. 
(Reference  temperature  25  °C,  and  pressure  at  1  bar 

Fillers  and  Tschoegl,  1977) 
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These  three  rubbery  materials  provide  numerical  estimates  of  how  much  the  time  scale 
changes  as  a  result  of  pressure  (in  kbars).  This  is  summarized  in  Table  II 


P  =  2  kbar 

P  =  4  kbar 

Hypalon  40 

2.5 

8 

Viton  B 

7 

- 

Neoprene  WB 

5 

11 

Table  II:  Shift  factors  for  pressure  increases  for  three  elastomers 


It  is  apparent  that  very  significant  stiffening  can  result  in  a  high  pressure  environment.  To 
the  best  of  our  knowledge,  this  response  is  valid  under  isostatic  as  well  as  dynamically 
changing  conditions.  In  the  latter  case  it  would  be  necessary,  in  principle,  to  prescribe  a 
specific  pressure  history  and  execute  the  integration  in  equation  (31-33)  to  evaluate 
elastomer  “stiffness”  in  connection  with  the  constitutive  law  outlined  for  small  strains  in 
equation  (23-24  and  sequel).  Just  to  keep  matters  in  perspective,  recall  that  in  the  material 
examined  in  the  present  work,  that  a  shift  of  8  decades  to  longer  times  (stiffening)  arises 
from  a  lowering  of  the  temperature  by  about  50  °C,  which  brings  the  material  close  to,  if 
not  directly  to,  the  glass  transition  temperature.  Thus  under  high  pressures  the  material 
responds  highly  viscously  or  in  a  leathery  manner. 


2. 2. 8. 2  Need  for  data  for  subsequent  analyses. 

To  date  the  above  referenced  data  are  virtually  the  only  ones  available,  or  certainly  the 
most  complete.  More  recent  work  by  Qvale  and  Ravi-Chandar  (2004,  see  section  1.1.6)  is 
more  limited,  but  points  in  the  same  direction. 

What  is  most  important  in  the  present  characterization  effort  is  to  recognize  that  both  the 
pressurization  process  as  well  as  the  high  rate  (short  time  scale)  leads  additively  to 
material  “stiffening”  which  corresponds  to  a  time  scale  that  is  not  covered  by  the 
measurements  in  this  report. 


3.  APPLICATION  TO  ARBITRARY  TEMPERATURES  AND 
DEFORMATION  HISTORIES 

The  most  immediate  question  relates  to  the  interpretation  of  the  relaxation  modulus  in 
figure  17  when  a  use  temperature  different  from  the  reference  temperature  of  ~0  °C  us  of 
interest.  Assuming  that  the  application  of  this  material  is  at  15  °C,  the  master  curve  of 
figure  17  has  been  replotted  in  figure  28;  this  amounted  to  shifting  the  curve  to  shorter 
log-times  by  an  amount  of  1.074  logarithmic  units  which  corresponds  to  the  time 


40 


Log(E*Tref?T)  MPa 


3.0 


2.8 

2.6 

2.4 

2.2 

2.0 

1.8 

-12  -11  -10  -9  -8  -7  -6  -5  -4  -3  -2  -1  0  1  2  3  4  5  6 

Log  t/a,  seconds 

Figure  28:  Relaxation  modulus  at  15  °C.  (The  ordinate  adjustment  from  figure  17 
is  only  5%;  on  a  log  scale  this  means  a  change  of  2%and  has  not  been  added) 
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temperature  shift  computed  from  the  WLF  equation  (  ).  Because  the  ratio  of  the  absolute 
temperature  for  15  °C  and  the  reference  temperature  is  278/273  is  small,  the  adjustment 
of  the  data  along  the  (logarithmic)  ordinate  is  very  small  and  has  not  been  incorporated 
into  that  figure. 


3. 1  VARIABLE  TEMPERA  TURE  AND  PRESSURE 

We  recall  the  three-dimensional  equations  of  (linear)  viscoelasticity 

<9e,(x,£) 

,S|,(x,0  =  fi«'-O  dt 

0 


(35) 


0 


cje(x,<*)-aAT (*,<*)] 


as  well  as  the  definitions  of  the  temperature/pressure  reduced  time(s) 


or 


with  the  shift  factor  aT,P  is  given  by  equations  (33)  and  (34). 


(36) 


(37) 


(38) 

(39) 


3.2  APPROXIMA  TIONS  FOR  FINITE  DEFORMA  TIONS. 

While  there  are  numerous  formulation  issues  associated  with  large  deformations,  the 
primary  question  in  the  present  context  arises  with  regard  to  interpretation  of  the 
nonlinear  behavior  in  the  uniaxial  data.  While  an  initial  approach  would  be  to  follow  the 
free  volume  concept  offered  by  Knauss  and  Emri  (1981,  1987),  or  the  stress-induced 
nonlinear  approach  discussed  by  Schapery  (1968,  1969),  these  concepts  reduce,  in  the 
limit  of  near-rubbery  behavior  (long-time  response)  to  a  multiplicative  decomposition 
into  the  viscoelastic  component  and  a  nonlinear  strain  function  as  long  as  the  defonnation 
follows  a  constant  rate  history. 
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Discussing  this  in  terms  of  uniaxial  behavior  would  yield  with  Sn=R-t,  as  a  simple 
approximation 

an{t)  =  Rsn)R\E{t-^  (40) 

o 

where  f  (sv] )  is  an  experimentally  detennined  function.  This  relation  states  simply  that 

the  time  dependencies  summarized  in  the  viscoelastic  integral,  and  the  strain  dependence 
in  the  function/  If  one  defines  the  viscoelastic  secant  modulus  or  constant  rate  modulus 
ER(t)  as 


E„(t)  =  -\E(l^£)d£ 
t  0 

then  (41)  may  be  written,  upon  noting  that  Sx  x  =  Rt,  as 


(41) 


o-n(t)  =  f(£ n H>n  •£«(»]  (42) 

This  relation  may  be  used  to  correlate  constant  strain  rate  (constant  deformation  rate)  data 
as  obtained  from  Instron  or  MTS  instruments.  Note  that  the  quantity  in  brackets 
represents  the  linearly  viscoelastic  response  in  a  uniaxial,  constant  strain  rate  test,  such  as 
would  be  computed  from  the  relaxation  data  in  figure  17  or  28. 


3.3  APPROXIMA  TIONS  FOR  MONOTONICALL  Y  CHANGING 
DEFORMATIONS 

The  following  is  a  simple  method  for  estimating  stress  and  deformations  in  complex 
situations.  It  gives  erroneous  results  if  cyclical  deformations  are  involved,  or  if  rates 
change  sign.  The  crux  of  the  method  is  to  detennine  a  reasonably  good  estimate  for  the 
material  stiffness.  This  value  is  governed  by  the  temperature  (history),  pressure  (history) 
and  deformation  rate.  If  deformations  remain  small,  in  if  rates  are  monotonic,  a  decent 
approximation  for  the  stiffness  of  the  material  is  obtained  by  the  following  process: 

a)  develop  the  relaxation  master  curve  for  the  appropriate  temperature  (and 
pressure); 

b)  determine  the  time  scale  of  interest;  if  the  defonnation  is  continuous  choose  the 

final  time;  if  relaxation  is  involved,  choose  the  final  time. 

c)  Determine  the  modulus  for  the  time  chosen;  this  value  is  a  reasonable  first 
estimate  for  stress  and  deformation  analysis. 

The  reason  that  this  method  is  useful  is  that  the  modulus  changes  significantly  only  if 
deformation  rates  change  by  an  order  of  magnitude:  it  is  a  modulus  that  is  relatively 
insensitive  to  details  of  the  deformation  history  (which  is  detennined  by  way  of  an 
integral) 

The  results  of  such  a  computation  renders  the  right  order  of  magnitude  of  quantities,  and 
is  probably  within  +/-15  %  of  the  desired  linearly  viscoelastic  result. 
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